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Abstract 

The multiplicity distributions in high energy e + e~ annihilation are described within 
the Two Stage Model. It is shown that oscillations in sign of the ratio of facto- 
rial cumulant moments over factorial moments of increasing order are confirmed. 
Parameters of parton stage and hadronization of this model are analyzed. 
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1 Introduction 



Multiparticle production (MP) is one of important regions in high energy physics pQ. 
Modern accelerators had made possible the intensive and detailed study of multiparticle 
processes. Quantum chromodynamics (QCD), developing theory of high energy physics 
[2] and a lot of phenomenologic models are tested by the process of MP. We should use 
phenomenologic models on account of absent total theory of QCD. 

MP is begun at high energy. Among of all producing particles we can observe a lot of 
hadrons. So, number of producing particles can be 50 and more in e + e~ annihilation [Hj. 

On one hand we want to know high energy physics, but on the other hand the increase of 
the inelastic channels is making difficulty in the description of this process with customary 
methods by using of the diagram technics. We have similarly situation as we had at the 
time of developing the thermodynamics and statistical physics. The total analysis of MP 
process is very difficult. The consideration of such reactions begins from the analysis of 
behavior of charge multiplicity. 

Definition: Multplicity is the number of secondaries n in process MP: 

A + B -> a x + a 2 + . . . + a n (1) 
Definition: Multiplicity distribution (MD) P n is the ratio cross-sections a n to a = On 

n 

P n = cr n /cr. 

This quantity has such sense: the probability of producing of n charge particles in process 
©• 

The description MP with using statistical methods is based on a lot of particles in 
this process. We can construct common quantities like mean values, moments of MD, 
correlations and so on. 

The history of study MP is very interesting. It connects on one hand with the developing 
theory and on the other hand with increase energy of accelerators. We have stopped on 
the studing MP in e + e~ - annihilation at high energy. At accordance with theory of strong 
interactions QCD this MP can be realize through the production of 7- or Z°-boson into 
two quarks: 

e + e - (£0/ 7 ) _ m (2) 

We can analyze the general features of hadronic decays up to the highest available 
center-of-mass (cm.) energies. In the last years LEP produced e + e~ -collisions of cm. 
energies at 172, 183 and 189 GeV. The charge particle multiplicity can be more than 60 

m 

QCD can describe process of fission partons (quarks and gluons) at high energy, because 
the strong coupling a s is small at that energy. This stage can be called as the stage of 
cascade. After fission, when partons have not high energy, they must change into hadrons, 
which we can observe. On this stage we shouldn't apply QCD, because a s is not small. 
Phenomenological models are used for description of hadronization (transformation partons 
(quarks and gluons) into hadrons). 
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It is expected on the stage quark-gluon cascade perturbative QCD will be applied 
jUE]. Certain features of the predictions at the parton level are expected to be insensitive 
to details of the hadronization mechanism. They can be tested directly by using hadron 
distributions 

The hadronization models are more phenomenological and are built on the experience 
gained from the study of low-p T hadron collisions. It is proposed to use Two Stage 
Model (TSM) [H] for description MD and other characteristics of e + e~ annihilation. It is 
usually considered that the producing of hadrons from partons is universal one. We can 
use these MD for calculation factorial and cumulant factorial moments and oscillations of 
their ratio and so on. 

Investigation of MP led to discovery of jets. Jets phenomena can be studied in all 
processes, where energetic partons are produced. The most common ones are in e + e~ an- 
nihilation, deep inelastic scattering of e, fi or v on nucleons and hadron- hadrons scattering, 
involving high-px particles in final state. 

The e + e~-reaction is simple for analysis, as the produced state is pure qq. It is usually 
difficult to determine the quark species on event-by-event basis. The experimental results 
are averaged over the quark type. Because of confinement the produced quark and gluons 
fragment into jets of observable hadrons. 

2 Two Stage Model 

Parton spectra in QCD quark and gluon jets have been studied by Konishi K., Ukawa A. 
and Veneciano G. 4J. Working at the leading logarithm approximation and avoiding IR 
divergences by considering finite x, the probabilistic nature of the problem has been estab- 
lished H]. 

At the study of MP at high energy we decided to use idea A. Giovannini [7j for describing 
QCD jets as Markov branching process. Giovannini proposed to interpret the natural QCD 
evolution parameter as the thickness of the jets and their development through subnuclear 
matter as Markov process. The Markov nature of elementary process ("gluon fission", 
"quark bremsstrahlung" and "quark-pair creation") leads to unsuspected properties of 
QCD jets. The stochastic approach gives clean complete solutions for the parton MD. 
QCD evolution parameter is 

Y= 2^b ln[1 + abln( ^ )] ' (3) 

where 2nb = |(lliVc — 2Nf) for a theory with Nc colours and Nf flavours. 

Three elementary process contribute to the quark or gluon distributions into QCD jets: 

(1) gluon fission; 

(2) quark bremsstrahlung; 

(3) quark pair creation. 

Let AAV be the probability that gluon in the infinitesimal interval Ay will convert 
into two gluons, AAY be the probability that quark will radiate a gluon, and BAY be the 
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probability that a quark-antiquark pair will be created from a gluon. A, A, B are assumed 
to be Y-independent constants and each individual parton acts independently from others, 
always with the same infinitesimal probability. 

Definite the probability that m g gluons and m q quarks will be transformed into n g glu- 
ons and n q quarks over a jet of thickness Y and call it P mg ,m q ;n g ,n q (Y). The probability 
generating function for a gluon jet and quark jet will be, respectivily 

oo 

G(u g ,u q ;Y) = £ P 1>0 . ng!nq (YX°u n q «, (4) 

n g ,n q =0 
oo 

Q(u g ,u q ; Y) = £ P 0A . ng;nq (Y)u n g °u n q «. (5) 

n g ,n q =0 

Due to the independent action of the individual partons it can be shown through straight- 
forward that 

oo 

E Pm g ,m q ;n g ,n q (YX"u n q " = [G( Ug , U q ;Y)] m » [Q(u g , U q] Y)] m « . (6) 
n g ,n q =0 

Eq. (JHj) says that from a probabilistic point of view the total parton population (m g gluons 
and m q quarks) evolving through thickness Y behaves as (m g + m q ) independent parton 
population, each with one initial quark or gluon. Since, the process is homogeneous in Y 
the transition probabilities obey Chapman- Kolmogorov equations: 

oo 

Pm g ,m q ;n g ,n q {Y ~\- Y ) — Pm g ,m q ;l g ,l q {Y)Pl g ,l q ;n g ,n q {Y ) • (^) 

lg,l q =0 

For the gluon jet one obtains 

oo 

Po,l;n g ,n q {Y + Y') = Pl,0;l g ,l q (^)-Pi SJ / B ;n s ,n g (Y') (8) 

l g ,l q =0 

and for a quark jet 

oo 

Po,l;n g ,n q (Y + Y') = P 0)1 . lg) i q (Y)Pi g} i q . ngtnq (Y'). (9) 

lg,l q =0 

Then (see eq.®,©) 

G(u 9 , u q - Y + Y') = G[G(u g , u q - Y'), Q(u g , u q - Y')- Y], (10) 

Q(u g , u q - Y + Y') = Q[G(u g , u q , Y'), Q(u g , u q - Y')- Y]. (11) 
Making the substitution Y' — > AY, we obtain 

dG(u a ,u a ;Y) dG (n s, . dG (n s, . . 

ay = ^ ro< "("""») + a^"^^ < 12 » 



3Q(u g ,u q ;Y) dQ () dQ () 

BY = dC^""*) + dZ W {U - U ^ (13) 



"g 

where the infinitesimal generating function for quark and gluon jets 
(w^ 9 \ug,u q ) and w^ q \u g , u q ), respectively) are 



w^(u 9 ,u q ) = £ a^l m u^u^ = (-A-B)u g + Au 2 g + Bu 2 q , (14) 



m 9 ,m 9 =0 

Eqs. (jl2j) and (fT3j) recognize the forward Kolmogorov equations for the generation 
functions of the transition probabilities P mg , mq]ng ,n q (Y). The corresponding backward Kol- 
mogorov equations follow from eqs. (flOll 1|) 

BG 

— = [G(u g , u qi Y), Q(u g , u q ] Y)], (16) 

|| = w^[G(u g , u q , Y), Q(u 9 , u q , Y)]. (17) 
Recalling eqs. (fH^ and eqs. (fTB^) and (|T7jl become 

BG 

— = -AG + AG 2 -BG + BQ 2 } (18) 

^ = -AQ + AQG. (19) 

These eqs. were obtained in [7j. We can find the probability for quark (or gluon) to 
produce in the interval (Y + AY) n g gluons and n q quarks through process (l)-(3). It 
follows for gluon jet 

Pi,o;n g ,n q (Y + AY) = [1 - An q AY - An q AY - Bn g AY]P 1}0 . ng , nq (Y) 

+An q AYP 1Ang ^ nq + A(n g - l)AYP h0 . ng ^ nq (Y) 

+ B(n g + 1) AFP li0 . n9+ling _ 2 (F) + o(AY) (20) 
Letting AY — > we get the system of differential equations: 

dPi, -n^ q (Y) = ^ _ B ^ _ An ^p i Qng nt] 
+An g P 1 

,0;n g — l,n q 

(Y) + A(n g - 1)A 

,0;n 9 — l,n q (Y) 

+ B(n g + l)P 1 , -n g+ l,n q -2(Y). (21) 
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Finding the explicit solutions both in terms of generating functions (|18| h ()19 )) or of 
the exclusive cross-sections (|20|) is not easy. Approximate solutions can be obtained for 
particular case B = 0, A ^ A ^ 0. 

The meaning of this approximation is next: we don't allow gluons to split into quark- 
antiquark pair. Eqs. for gluon and quark jets in this case are 



dPi,o-,n g ,o{Y) 



dY 

dPoA;n g ,l\Y) 

dY 



A(n g - l)Pi, O ;n e -i,o0O - An g P l70 . ngfi (Y), 



(22) 



-AP 0tl;ngtl (Y) - An g P 0)1 . ngA (Y) 



+ + An g P 0tl . nail (Y) 



(23) 



with initial conditions P ,i;0,i(°) = h A,o;i,o(0) = 1, Pi,o;n g ,o{ Q ) = and ^b,i;n ff ,i(°) = 
Vn 9 > 1. 

For gluon jet we have 



n„-l 



Pl,0;n 3 fi{Y) 



<n g > 



< n g > t 



(24) 



with a average gluon multiplicity < n g >= e AY . For quark jet ()23j) one obtains MD 



P ,i;o,i(Y) = e 



-AY 



Po^(Y) = ^ + 1 )-Cy + ^-D e -ar (1 _ e -Ayy^ 



n g \ 



(25) 
(26) 



where \x 



The average gluon multiplicity will be < n g >= fi(e AY — 1) and the 



normalized exclusive cross section for producing n g gluons 

r(0,9) 



Otot 

/i(/i + 1) ...(/! + n g — l) 



Po,l;n g ,l[Y) 



n g \ 



< n y> 
< n g > +fi 



< n g > +/i 



The generating function will be given by 

Q X U n g °U q P ,l;n g ,l{Y) 



u„ 



n o =0 



-AY 



l-u g {l-e-*v) 



(27) 



(2f 



Eq.(|2Tj) is Polya-Egenberger distribution, where \i is not integer. 

In Two Stage Model jHj we took (|25|l for description of parton stage and added for final 
stage hadronization supernarrow distribution of Bernulli (binomial distribution). We chose 
it based oneself on analysis experimental data in e + e~- annihilation at 9 - 22 GeV. Second 
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correlation moments were negative at this energy and it was one, which could describe 
experiment. 

We consider that hypothesis of soft bleachment is right. We added stage of hadroniza- 
tion to parton stage with aid of their factorization. MD in this process can be writing 



(29) 



where P^ - MD for partons (|27|) . P^(m, s) - MD for hadrons, are produced from m partons 
on the stage of hadronization. 

In accordance with TSM the stage of hard fission of partons is described by negative 
binomial distribution (NBD) for quark jet 



k p (k p + 1) . . . (k p + m — I) 
ml 



m \ 



km 



m + k v 



k p + m / 



(30) 



where k p = A/ A, m = J2 m Pm- We are neglecting process (3) quark pair production 

m 

(B = 0). Two quarks fragment to parton independently each other. Total MD of two 
quarks is equal to ([50]) too. MD P^ and generating function for MD Q p (s, z) are 



1 d n 



ml dz r 



Q p (s,z) 



Qm( s > z ) 



l + r (l-z) 

tvp 



2 = ' 
"I &v 



From TSM MD of hadron are formed from parton are described in form 



P. 



H 



n 



N v -n 



p \N, 



N„ 



(31) 
(32) 

(33) 



(C^ - binomial coefficient) with generating function 



Qv 



1V 



N v 



(34) 



where and N p (p = q, g) have sense of average multiplicity and maximum secondaries 
of hadrons are formed from parton on the stage of hadronization. MD of hadrons in e + e~ 
annihilation are determined by convolution of two stages 



x P d r ' 



Pn(s) = E P, 



dz 1 



{Q 



H\2+m\ 



2=0 



(35) 



Further we assume next simplification for second stage: ^- « jf- . Let a = jf-, N = N t 



n • = n q . Then 



«f = (i+7r«*-i))° ■ 

From(j2SJ) we have 

Pn(s)= E P ^ C 2 + arn)N [j^J i 1 ~ n) " ( 36 ) 

For comparing with experimental data the normalized factor Q was introduced into (}36j) 
and a number of summarized gluons were restricted by M g - maximal number of gluons 
are giving cascade 

M g /rT h \ n ( w h\( 2 + am ) N - n 

Pn(s) = nJ2prc? 2+am)N l-j (i--J . (37) 

The results of comparison of model expression MD (f3Tj) with experimental data [E| 
represented in the Table 1. Apparently that behavior of parameters of TSM is giving the 
reasonable picture of e + e~ annihilation process. The average multiplicity of gluons m g 
formed on fission stage changes from 3 - 4 at 22 GeV to 15 - 16 at 189 Gev. This increase 
takes place with the growing energy, but at the some energy this value accepts a big one 
(at 43.6 GeV). 

From QCD is following that the parameter k p must decrease at more high energies. 
The parameter of parton spectra k p changes: from 140 at low energy (22 GeV) up 5 - 6 at 
172 - 189 GeV. QCD are giving limiting value k p is equal to 1. 

The most interesting dependencies are discovered for parameters of hadronization n^, 
N q and a . The parameter N q determines maximum number of hadrons, which can be 
formed from quark on stage of hadronization. It changes unusually, it can be 5 - 6 or 8 - 
15 at different energies, so it equal about 9 at 172 GeV, but equal 54.6 at 183 and again 11 
at 189 GeV. After that analysis we can say about the existing some special energies when 
process hadronization changes (at those energies value of x 2 is changed unconsiderably) . 

At the same time the parameter of hadronization (mean number of hadrons formed 
from quark on stage of hadronization) discovers insignificant growth with increasing of 
energy. It is changed from 3 - 4 at 22 GeV up 4 - 5 at 172 - 189 GeV. The existing models 
give quantity rig equal two (it is observed really at energies lower than 40 GeV). Such 
behavior of parameter may be connected with the growth of spectrum of mass hadron 
states (appearance of new mass states with increase of energy). It is interesting that the 

ratio jf- (the parameter of PBD) is approximately regularity and is equal to 1/2. It should 
mark breaking such behavior at some energies (43.6, 183 GeV) again. 

Parameter a is stayed approximately constant and equal to 0.2 in order case and become 
smaller than 0.1 at special cases (43.6 and 183 GeV). The number of hadrons are produced 
from quark jet on the stage of hadronization is more than from gluon jet. The normalized 
factor Q is stayed constant and equal to 2. 



7 



Figure 1: Multiplicity distributions P n and H q as a function of q at Z°. 



3 Oscillation of moments in MD 

At the recent years it was shown jH] that the ratio of factorial cumulative moments over 
factorial moments changes sign as a function of order. We can use MD formed in TSM for 
explanation of this phenomena. 

The factorial moments can be obtained from MD P n through the relation 

oo 

F q = J2 n{n - 1) . . . (n - q + l)P n , (38) 

n=q 

and factorial cumulative moments are found from expression 

K^F.-^C^K^F,. (39) 
i=i 

The ratio of their quantities is 

H q = K q /F q . (40) 
We can use the generating function for MD of hadrons (j3*U|) in e + e~ annihilation G(z) 

G{z) = E Pi[Q?(zTQ 2 q (z) = 

m=0 



Q9(Q»(z)rQ 2 Jz). (41) 



We are calculating F q and K q in TSM using ()41|) 

1 d°G 



K n 



q n*(s) dzi 
1 d g \nG 



2=0 



' q W(s) dz* 
The expression f!41|) for G(z) after taking a logarithm 

lnG(s, z) - 



(42) 
(43) 



m 



z=0 



fc p ln[l + r (l-Q£)]+21nQ 

Kp 



and the expansion to series in power on will be 



lnG(s,z) = k P Y^ 



rn 



q: 



\ \m + k p 



rn 



+ 21nQ. 



(44) 



Inserting Q g into (jUJ) 



In G{s, z) = k p ^2 [ — 



m 



m=0 



m + k v 



1 

m 



1 + Tf( z ~ ' 



amN 



+ 



n 



+2tfln[l + —(*-!)], 



and using (|4*3j) we obtain 



V m=l 



1> 



rn 



[am 



N 



m + k v 



1 

m 



1; 



.(g-i)!' 



n 



-i 



n is 



where n(s) is the average multiplicity hadrons in process ( 
(PJ) 

F g = ^ (2 + am) (2 + am - h . . . (2 + am - 

m=0 ^ V ^ V 

with P m equal (|50|) . 

The sought-for expression for if g will be 



(45) 

It is possible to find F 9 using 



n 



n(s) 



(46) 



E M m ( am -£)... (am - 2^)(=^ ; - n -x ; 



— 2f— 11* 



(9-1)! 



m=l 



E (2 + am) (2 + am - £) ... (2 + am - ^)P, 

m=0 



(47) 



where fix is the normalized factor. The comparison with experimental data ^U] shows that 
(}4Tj) is describing the ratio of factorial moments. It is seen minimum at q=5. 

The immediate calculations H q based on (J38|) - f4*0|) with using MD (j37j) are giving very 
good description of the oscillation value of H q (x 2 rs 2). Significant oscillations begin 
near region producing of Z° and can be explain by not integer values of parameters of 
hadronization N q and N g = aN q or by convolution of wide NBD and narrow PBD. 
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4 Conclusions 



It is shown that TSM does not contradict to the experimental data on MD and the oscil- 
lations ratio of factorial moments (small x 2 ) (Figure 1). 



Table 1. Parameters of TSM. 



y/s GcV 


m 


kp 


N 


n h 


a 


n 


14 


.08 


2.4* 10 8 


27.7 


2.87 


.97 


2 


22 


3.01 


4.91 


20.2 


4.34 


.2 


2 


34.8 


6.58 


6.96 


12.5 


4.1 


.195 


2 


43.6 


10.9 


39.4 


6.1 


2.43 


.386 


2 


91.4 


10.9 


7.86 


11.2 


4.8 


.226 


2 


172 


20.1 


9.11 


9.17 


4.34 


.195 


1.98 


183 


13.2 


1.48 


54.6 


8.9 


.086 


2.06 


189 


15.1 


6.9 


11.6 


5.15 


.215 


2.01 
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